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Abstract: This paper presents an acoustic noise cancelling technique using an inverse 
kepstrum system as an innovations-based whitening application for an adaptive finite 
impulse response (FIR) filter in beamforming structure. The inverse kepstrum method uses 
an innovations-whitened form from one acoustic path transfer function between a reference 
microphone sensor and a noise source so that the rear-end reference signal will then be a 
whitened sequence to a cascaded adaptive FIR filter in the beamforming structure. By 
using an inverse kepstrum filter as a whitening filter with the use of a delay filter, the 
cascaded adaptive FIR filter estimates only the numerator of the polynomial part from the 
ratio of overall combined transfer functions. The test results have shown that the adaptive 
FIR filter is more effective in beamforming structure than an adaptive noise cancelling 
(ANC) structure in terms of signal distortion in the desired signal and noise reduction in 
noise with nonminimum phase components. In addition, the inverse kepstrum method 
shows almost the same convergence level in estimate of noise statistics with the use of a 
smaller amount of adaptive FIR filter weights than the kepstrum method, hence it could 
provide better computational simplicity in processing. Furthermore, the rear-end inverse 
kepstrum method in beamforming structure has shown less signal distortion in the desired 
signal than the front-end kepstrum method and the front-end inverse kepstrum method in 
beamforming structure. 

Keywords: innovations; whitening; ANC; beamforming; inverse kepstrum; adaptive FIR 
filter; system identification; acoustic transfer function 
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1. Introduction 

During the last five decades, noise cancelling and signal enhancing techniques have been developed. 
The techniques are fundamentally based on spectral subtraction, cepstrum and complex 
cepstrum methods using single-microphone sensor, and ANC and beamforming methods using 
multiple-microphone array sensors. 

Since research on echo cancellation using adaptive filters and two-microphone sensors started in 
1965, the adaptive filtering technique has been used as a solution tool for signal enhancing and noise 
cancelling schemes [1]. In 1975, Widrow et al. [2] proposed an FIR least mean square (LMS) 
algorithm-based ANC method using two microphones. This method requires that the primary 
microphone takes signal plus noise and the reference microphone takes noise alone, whereby the 
adaptive filter estimates noise statistics from the reference microphone and then uses them to minimize 
the output power in a minimizing mean square error (MMSE) calculation. From the theory, the 
application shows a practical problem due to the difficulty of separation between the period of the 
noise alone and the period of signal plus noise. This causes signal leakage into the reference 
microphone, which makes the adaptive filter estimate noise as well as the desired signal, hence results 
in limitations in the maximum cancellation from the output signal-to-noise ratio (SNR) with 
signal distortion. 

Pulsipher et al. [3] have analyzed an unknown system used for the identification of acoustic path 
transfer functions between two-microphone sensors and noise sources from the ANC method and 
investigated the adaptive filter estimates of the ratio of acoustic path transfer functions for the 
correlated noise and found out that the problem comes from the nonminimum phase component of the 
ratio of acoustic transfer functions. To resolve it, the application uses a large amount of adaptive 
filter weights. 

Harrison et al. [4] have introduced a new approach in the ANC method, whereby they use a noise 
estimating technique of a small separation between two-microphone sensors with the use of a voice 
activity detector (VAD) during the noise alone period. The result has shown that it could significantly 
reduce the amount of adaptive filter weights required for noise cancellation while minimizing the 
presence of reverberation. Nevertheless, the ANC method shows a limitation in the maximum 
cancellation in SNR because this maximum cancellation is related to the coherence level of the noise 
between two-microphone sensors [5], whereby it requires a significant coherence for even modest 
noise cancelling performance. The application of close and direct application of a desired signal [6-8] 
in front of two microphones with the use of adaptive filter using VAD during the noise alone period 
may reduce signal distortion. 

The beamforming technique has been introduced to maximize signal directivity, therefore it 
increases the performance in SNR. This method may require many microphones with accurate phase 
alignment among the microphones array and hence, a computational complexity in processing is 
expected. To increase the performance in SNR, the technique has been developed with the use of 
adaptive filter and also VAD [9]. 

The cepstrum processing technique [10] may provide a solution for signal separation, but with the 
practical limitation due to minimum phase information. The complex cepstrum method [11] may give 
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the solution for the minimum phase problem, where this minimum phase information can be directly 
estimated from the power spectrum. 

The kepstrum [12,13] is similar to the complex cepstrum due to the fact that the minimum phase 
spectral factor can be directly obtained from a power spectrum estimation and it is equivalent with 
logarithmic minimum phase transfer function. The kepstrum method has been used as a system 
identification technique for unknown systems from the acoustic path transfer functions between 
microphone sensors input and the noise source in beamforming structure [14,15], with a phase 
recovering technique from the causal kepstrum domain, where the front-end kepstrum method has 
produced an improved SNR performance in a different input SNR for real-time processing in 
reverberant room environments. It has also been tested in different locations on a noise source using 
nonstationary music sound, tuned to a radio station. With a kepstrum method, a technique for 
preventing signal distortion has been used with the modified application [6-8] to a desired signal, and 
also modified application [4] to an adaptive filter with the use of VAD for differentiating the periods 
between signal plus noise and noise alone. 

In addition, the kepstrum has a distinction in the case of signal plus noise, where the logarithmic 
minimum phase transfer function becomes the minimum phase kepstrum spectral factor and it can be 
represented as a Kolmogorov [16] power series expansion. Furthermore, the signal and noise may be 
implemented in an innovations-based form, where it was originally discovered by Kalman and 
Bucy [17] and it may then be applied to an infinite impulse response (IIR) Wiener filtering structure [18]. 
For the innovations-based approach, it may be represented as an output of normalized minimum phase 
spectral factor from innovations white noise input and has been used by Moir and Barrett [19]. 

By applying innovations-based whitening application in an ANC structure, it has been investigated 
in a simulation test, where it has been applied with the use of a FIR normalized least mean square 
(NLMS) algorithm for noise cancellation [20] and also with a FIR recursive least squares (RLS) 
algorithm [21], where it has been found that the application of innovations-based inverse kepstrum to 
cascaded adaptive filter gives a stable and causal system because all poles and zeros of the system are 
converted into unit circles due to the fact that the whitening application in the reference input works 
as all-pass filter so that it allows only one acoustic path transfer function to be considered as the 
unknown system. 

For the real-time processing using adaptive RLS filters, the innovations-based whitening 
application has been applied as a front-end application in beamforming structure with rear-end 
zero-model FIR RLS filters and it was found that it gives better noise cancelling performance than a 
pole -zero model IIR RLS filter in an ANC structure [22]. 

From the previous studies [14,20-22] based on modified applications [4,6-8], it has been found that 
there are important features to be further investigated to verify the performance. In this paper, by 
considering: (1) signal distortion in a desired signal on the instant application of noise statistics in 
reverberant environment, (2) noise reduction on noise characteristics with nonminimum phase 
components and its consistency on its inverted acoustic path transfer function, and (3) the use of small 
amount of adaptive FIR filter weights for the real-time processing, hence for the fast convergence in 
estimate on noise statistics, it is analyzed in both ANC and beamforming structures. The inverse 
kepstrum method is then applied to the rear-end (of sum-and-subtract function) in beamforming 
structure to verify the performance by comparing signal waveforms in the time domain, spectra in the 
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frequency domain as well as estimated coefficients arrays of inverse kepstrum and weights arrays of 
adaptive FIR filter and its pole-zero placements in noise statistics. Furthermore, the rear-end inverse 
kepstrum method is also compared with the front-end kepstrum method [14], which uses identification 
of two paths with the ratio of overall acoustic path transfer function, and also with front-end inverse 
kepstrum method using a whitening application [22]. 

2. Analysis of Innovations-Based Inverse Kepstrum 

This section describes the analysis of cepstrum and complex cepstrum (kepstrum) with the relation 
of minimum phase kepstrum, and also for its whitening application (inverse kepstrum). It shows that 
the minimum phase kepstrum coefficients may be obtained from the logarithm of the minimum phase 
transfer function (Section 2.1) and also from the logarithm of the minimum phase spectral factor 
(Section 2.2). For the signal and noise, it shows that it can be represented as an output of normalized 
minimum phase spectral factor from innovations white noise input (Section 2.3). Based on this, it 
shows that logarithm of inverse minimum phase transfer function may be implemented as the 
innovations form of the normalized minimum phase kepstrum spectral factor for the whitening 
application (Section 2.4). 

2.1. Analysis of the Minimum Phase Transfer Function 

It is known that the causal transfer function can be expressed by Schwarz's classical formula [23] 



as: 



dl , \z\ < 1 



(1) 



where X as the integration variable and z = re jw . 

This Equation gives the causal transfer function H + (z) whose real part on the unit circle is H R (w). 
Based on this, the phase information can be recovered by Hilbert's transform relation. The logarithm 
of minimum phase transfer function \ogH M (z) can be written as: 



\ogH M {z) = ^)\og\H M {X)\ 



1 + e JA z 



l-e M z 



dX, \z\ < 1 
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where log \H M (X)\ is the magnitude part of minimum phase logarithmic transfer function. 

This indicates that the minimum phase transfer function may be expressed in terms of Schwarz's 
formula and hence the minimum phase information can be recovered from Hilbert's transform relation. 

By defining that z = e JW , it can be described as magnitude and phase term: 
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From Equation (4), the magnitude of logarithmic minimum phase transfer function is: 
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and the phase of the logarithmic minimum phase transfer function is: 
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This shows that Equation (5) shows an even cepstrum function, hence minimum phase kepstrum 
coefficients can be processed in the cepstrum domain by multiplying by two (2) in the positive time 
series, except for the first zeroth coefficient in the time series. 

2.2. Analysis of Minimum Phase Spectral Factor 

From the power spectral density O(z), it can be represented as causal spectral factor H + (z) and 
anticausal counterpart H~(z) as: 

cD(z) = H + (z)H(z) = H\z)H + (z 1 ) (7) 
Let z = z" 1 , then it follows that 0(z _1 ) = H + (z~ l )H + (z) 

It follows that 3>(z) = cD(z -1 ). By defining Z = e JW , we now have a logarithmic power spectrum 



as: 
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As a result of the symmetry property of k n , it can be expressed as: 
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(12) 



Furthermore, since are real, only a half portion in length can be considered as: 
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From Equation (13), it shows that the kepstrum coefficients can be processed in the causal kepstrum 
domain by halving the first zeroth coefficient with the remaining coefficients truncated in size 
to(iV72)-l. By truncating in size less than(A^/2)-l , kepstrum now becomes complex cepstrum, 
which is an approximation of the theoretical mathematical construct [19]. 
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2.3. Analysis of Minimum Phase Kepstrum Spectral Factor on Signal and Noise 

In the case of random signal plus noise, it can also be represented as innovations-based form. From 
Figure 1 (upper part), it shows an equivalent relation of output from the inputs between white noise and 
innovations white noise as Equation (14): 

® ** 00 = 00 + <D ss ( z ) = N M (z)N M (z-> )al + S M (z)S M (z~ l )a] = K {z)H~ n (z)af ( { 4) 

where o 2 v , a\ and cr 2 are the variances of the additive white noise, white noise input and white noise 
innovations process, respectively. N M (z) and S M (z) are coloured minimum phase transfer functions 
and N M {z~ v ) and S M (z~')are maximum phase counterparts. H + n (z) is the normalized minimum phase 
spectral factor, which has all its zeros inside \z\ = 1 whilst H~{z) is the counterpart, which has its zeros 
outside \z\ = 1 . 

Figure 1. Equivalence of outputs based on inputs of: (a) white noise model and 
(b) innovations white noise model [upper part: signal plus noise, lower part: noise alone]. 
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In the case of signal plus noise, the logarithm of each positive- and negative-sided transfer function 
becomes the kepstrum spectral factors of the z- transform spectral density and these are represented as 
a power series expansion. For the innovations-based inverse kepstrum approach, signal plus noise are 
represented as an output of normalized minimum phase spectral factor from the innovations white 
noise input. It may be applied to an optimum IIR Wiener filtering structure, where it has been defined 
by Kailath [18] as combination of two cascaded filters, a front-end whitening filter to generate the 
white innovations process and cascaded spectral shaping filter to provide spectral shaping function for 
the input signal. 

2.4. Analysis of Innovations-Based Inverse Kepstrum and its Application to Noise Signal Only 

For the application of noise alone, it can be estimated during the absence of desired signal 
(Figure 1 (lower part)). Therefore, the additive noise now can be represented as innovations-based form, 
as shown in Figure 2. 

Figure 2. Representation of noise as innovations-based whitening form. 
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From the minimum phase kepstrum spectral factor of Equation (15): 

K + (z) = \ogH + (z) =k 0 +k lZ - l +k 2 z~ 2 + (15) 

its exponentiation becomes a causal spectral factor as: 

H + (z) = Qxp[K + (z)] = Qxp(k 0 + k l z~ i +k 2 z~ 2 + ...) (16) 

By using the fact that minimum phase spectral factor allows its invertibility, it can be expressed as: 
[H + {z)Y l = exp[-JT (z)] = exp[-(£ 0 + k x z~' + k^ 1 + ...)] (17) 

This shows that the inverse of the minimum phase spectral factor can be obtained from 
the kepstrum exponential by multiplying by minus one (-1). It can also be represented in a 
normalized form: 

H + (z) = exp[k Q ]exp[K + n (z)] = exp[k 0 ]H + n (z) (18) 
logH + (z)=k 0 +logH + n (z) (19) 

In the case of the application of additive noise alone, l„ becomes V n . From the fact that 
x n =H M {z)v n and x n =H + n (z)i n , the relationship between minimum phase transfer function H M {z) and 
normalized minimum phase spectral factor H + n (z) is described as: 

H M {z) = H + n {z) i„ = H + n {z)v n = H + n (z)o c (20) 

where V n , i n and a e are white noise, innovations process and standard deviation of normalized 

innovations sequence respectively. 

By taking logarithm of Equation (20): 

log//„(z) = log/C(z) + logo; (21) 

From Equations (19) and (21) and with the fact that log H + (z) = \ogH M (z) : 

K=^gcr £ (22) 

a £ = exp [k 0 ] and a] = exp[2k 0 ] 

where o £ and a £ are standard deviation and variance of normalized innovations sequence 
respectively. 

A normalized innovations-based inverse kepstrum is then represented from logarithmic minimum 
phase transfer function or logarithmic minimum phase spectral factor, which described as: 

log [H M (z)]- 1 = log [H + n {z)Y l - \oga c = -K + n (z) - log u e (24) 

Therefore, it shows that logarithm of inverse minimum phase transfer function may be implemented as 
innovations form of normalized minimum phase kepstrum spectral factor (log [H M {z)y l = -K + n (z), 
assuming that standard deviation <j e = 1 ). 
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3. Inverse Kepstrum Processing and Adaptive FIR NLMS Algorithm 

For the whitening application, the inverse kepstrum is processed from the reference microphone x n . 
Therefore only the estimate from a single acoustic path transfer function is required. This may be 
compared with kepstrum processing, where the estimate of two acoustic path transfer functions from 
both the primary microphone d„ and reference microphone x n is required. From each input microphone 
sensor, periodograms are obtained from Hanning-windowed fast Fourier transforms (FFTs) from the 
two-microphone inputs as shown in Figure 3. 

Figure 3. Periodogram estimate from each input microphone sensors. 
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As a discrete estimate of the continuous power spectral density, it uses a modified weighted 
overlapped segment averaging (WOSA) algorithm and the auto-periodograms are processed from 50% 
overlapping Hanning-windowed FFTs in 2,048 frame size by using smoothing method (25), with the 
forgetting factor^ = 0.8 [20]: 

<M0 = FS>*,(i - 1) + (1 - P)X d (i)X' d (i) <!>„(i) = fiO„(i - 1) + (1 - 0)X x (i)X' x (j) (25) 

From (25), its logarithm needs to add Euler's constant to be unbiased due to bias in magnitude [24]. 
Hence, each kepstrum coefficient from the two-microphone inputs are found from the inverse FFT 
(IFFT) of the unbiased logarithmic auto-periodograms and then by subtracting two kepstrum 
coefficient vectors (k ln -k 2n ), we can get the kepstrum coefficients (k n ) for kepstrum processing. On 
the other hand, inverse kepstrum coefficients (k n ) are found by negating the kepstrum coefficients 
(-k 2n ) from the reference microphone. The processing difference between inverse kepstrum method 
and kepstrum method can be found in Figure 4. 

Figure 4. Block diagram for the comparison between inverse kepstrum processing ( - k 2n ) 
and kepstrum processing ( k ln - k 2n ). 
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This indicates that inverse kepstrum processing requires only a negative sign of the kepstrum 
coefficients, which can be obtained for the inverse of acoustic path transfer function from the reference 
microphone input Equation (26). The negated kepstrum coefficients are then normalized by dividing 
the zeroth kepstrum coefficient value: 
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log{l/# 2 (z)} 4*K\z) = -K 2 (z) 



(26) 



On the other hand, it is compared with the kepstrum method, which uses the ratio of the acoustic 
transfer function between two-microphone inputs as Equation (27). 

\og{H x {z)IH 2 {z)} ^ K(z) = K 1 (z)-K 2 (z) (27) 

Both the kepstrum and inverse kepstrum coefficients are transformed into the corresponding 
impulse response using recursive formula [12] and then it is convolved with reference input signal to 
get a refined new input signal as shown in Figure 5. 

Figure 5. The conversion procedure from inverse kepstrum ( k' n ) or kepstrum ( k n ) to 
impulse response (h n ), and its convolution with reference input signal (x„). 
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For the cascaded adaptive FIR filter, the NLMS algorithm [25] has been used and the weights are 
updated as: 



0.0001 + 1 \xj ( 28 ) 

where pi is step size (0 < /u < 2) , \\X n f is input power and the value of 0.0001 is used to prevent zero 
division. 



■ 1 2 . 



Table 1. Comparison of computational complexity in FLOPS on (i) inverse kepstrum 
processing, (ii) kepstrum processing and (iii) NLMS algorithm [20]. 



Algorithm 


Required processing 


FLOPS 


(i) Inverse kepstrum 


1 x WOSA Periodogram 


7An 0 g 2 (5.12/A/) 




4 x FFT / IFFT 


4(7V72)log 2 7V 




2 x Logarithm /Exponential 


2/V 1/3 (logA0 2 




Total Computation (*) 


0.057 x 10 6 


(ii) Kepstrum 


2 x WOSA Periodogram 


2tflog 2 (5.12/AO 




5 x FFT / IFFT 


5(7V72)log 2 7V 




3 x Logarithm /Exponential 


3N m (log N) 2 




Total Computation (*) 


0.08 x 10 6 


(iii) NLMS 


Real multiplication (A) 


3N 2 + 2N 




Iterations (B) 


20 x (A) 




Total Computation (**) 


(A) 0.12 x 10 6 (B)2.4 x 10 6 



Note: Total computation: (*) is based on N: 2048 frame size, (**) is based on 200 NLMS weights. 



For the comparison of processing among: (i) inverse kepstrum coefficients, (ii) kepstrum 
coefficients and (iii) NLMS algorithm based adaptive FIR filter weights, computational complexity in 
floating point operations per second (FLOPS) can be compared as shown in Table 1, where it shows 
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that inverse kepstrum processing gives the least computational complexity, and kepstrum processing 
shows less computational complexity than the ordinary adaptive FIRNLMS algorithm [20]. 

4. Inverse Kepstrum Method 

The inverse kepstrum method uses the whitening application from one acoustic path transfer 
function, H 2 (z) from the reference microphone during the noise alone period as shown in Figure 6(a), 
where the inverse kepstrum is to be analyzed as the front-end application to the adaptive FIR filter 
from the ANC structure, and the rear-end application to the sum-and-subtract function from the 
beamforming structure accordingly. It is also to be compared with the kepstrum method, which is 
based on identification of the ratio of acoustic path transfer functions, h { (z) and# 2 (z). During the 
signal and noise period, the noise estimate is applied to obtain the desired signal with a non-distortion. 
For the purpose, the method uses direct application of desired signal in front of two microphones as 
shown in Figure 6(b). Direct application of the desired signal can be found in the ANC structure [6] 
and also in beamforming structure [7,8]. 

Figure 6. (a) Acoustic path transfer functions as noise estimate and (b) application of 
desired signal to the two microphones. 
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4.1. Inverse Kepstrum Method as Front-End Application to ANC Structure 

In the ANC structure, identification as an unknown system is represented as the ratio of acoustic 
transfer functions, H(z) = H x (z)l H 2 (z)as> shown in Figure 7, where an ordinary adaptive FIR filter may 
be used to estimate the noise statistics during noise alone period. Its estimate is then applied to the 
signal and noise period to cancel the noise by using the estimated noise statistics with the condition 
that the desired signal should be only retained with no distortion. 

Figure 7. Analysis of identification of acoustic path transfer functions by an adaptive FIR 
filter in the ANC structure during the periods of (a) noise alone and (b) signal plus noise. 
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Secondly, the inverse kepstrum is applied in front of the adaptive FIR filter as shown in Figure 8(a), 
where the inverse kepstrum filter estimates a denominator part, \lH 2 (z) and the cascaded adaptive FIR 
filter estimates a numerator part, H l (z) from the ratio of overall transfer function. It is compared with 
the kepstrum method as shown in Figure 8(b), where the kepstrum filter estimates the minimum phase 
term only from the ratio of overall transfer function and the cascaded adaptive FIR filter estimates 
remaining term and hence it works as an all-pass filter. 

Figure 8. Application of (a) inverse kepstrum filter and (b) kepstrum filter to 
the ANC structure. 
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Assuming that each acoustic path transfer functions between two microphones and noise source are 
as given in (29), it is represented as the unknown system of (30). Its estimates are to be analyzed as 
shown in Figure 9, where the operation of the inverse kepstrum is to be compared with the kepstrum 
filter as a front-end application to the cascaded adaptive FIR filter. 

Figure 9. Block diagram for the operation comparison between an inverse kepstrum filter 
and a kepstrum filter as front-end application to an ANC structure. 
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As a simple example, we assume that one transfer function from acoustic path transfer functions is 
nonminimum phase term, such as: 

H x (z) = 1 + 2z~ x and H 2 (z) = 1 + 0.2z _1 (29) 

where H Y {z) is nonminimum phase transfer function. 

The unknown system is then described as the ratio of transfer functions and it is represented as: 
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H(z) = (1 + 2z _1 )/(l + 0.2Z" 1 ) = 1 + 1 - 0.36z~ 2 + 0.072z~ 3 - 0.0144z~ 4 (30) 

where it can be estimated by ordinary adaptive FIR filter. 

For the operation of the inverse kepstrum filter K (z) and adaptive FIR filter W x (z) , each one 

is estimated as: 

A"'(z) = l/(l + 0.2z _I ) W l (z) = l + 2z- ] (31) 



This indicates that the front-end inverse kepstrum K (z) estimates the denominator part and the 
cascaded adaptive FIR filter W x {z) estimates the numerator part from the ratio of overall transfer 
functions. On the other hand, for the operation of the kepstrum filter K(z) and adaptive FIR filter 
W 2 {z) , each one is estimated as: 

^(z) = (l + 0.5z _1 )/(l + 0.2z- 1 ) ^ 2 (z) = (l + 2z _1 )/(l + 0.5z- 1 ) (32) 

This indicates that the front-end kepstrum estimates the minimum phase term only from the ratio of 
overall transfer function, where nonminimum phase term is reflected to the minimum phase term by a 
reciprocal polynomial as z~"H(z~ l ) . The cascaded adaptive FIR filter works then as an all-pass filter. 

From (31) and (32), we now have found that both methods show the same result as (30). 
Now let us check with inverse of overall transfer function, such that: 

H.iz) = 1 + 0.2z _1 and H 2 (z) = l + 2z _1 (33) 

where H 2 (z) is nonminimum phase transfer function. 

This indicates that the transfer function now has a nonminimum phase term in the denominator 
polynomial from the ratio of overall transfer functions. 

The unknown system can be estimated from the ratio of transfer functions and it is represented as: 

H{z) = (1 + 0.2z-')/(l + 2z _1 ) = [(1 + 0.2z^)/(l + 0.5z^)][(l + O^z' 1 )/^ + 2z -1 )] 

= (1 + 0.2z"') /(l + 0.5Z" 1 ) = 1 - 0.3Z" 1 + 0. 15z~ 2 + (34) 

where it can also be estimated by adaptive FIR filter. The Equation (34) indicates that the 
nonminimum phase term of denominator part is converted to a minimum phase term for the operation 
of an adaptive FIR filter so that the cascaded remaining part of Equation (34) always works as an 
all-pass filter. 

The inverse kepstrum and cascaded adaptive FIR filter are estimated as: 

K\z) = l/(l + 0.5z~ l ) W l {z) = \ + 02z- 1 (35) 

and the kepstrum and cascaded adaptive FIR filter are estimated as: 

i:(z) = (l + 0.2z- 1 )/(l + 0.5z- 1 ) W 2 (z) = \ (36) 

This indicates that the estimates, Equations (35) and (36), from both methods are the same as 
Equation (34), which is the inverse of the overall transfer function. The nonminimum phase 
component from the polynomial numerator and denominator of the overall transfer function may 
frequently occur in reverberant environments [26] and causes a fluctuation in the spectrum due to the 
difference between Equations (30) and (34). To deal with this problem, it indicates that neither the 
kepstrum nor the inverse kepstrum methods provide a solution in the ANC structure. Alternatively, we 
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may solve this problem by swapping the position of the two microphones, but it is not a practical 
solution in a realistic environment. 

4.2. Inverse Kepstrum Method as a Rear-End Application in Beamforming Structure 

In a beamforming structure, identification as an unknown system is represented as the ratio of 
combined acoustic transfer functions, H(z) = 0.5(H i (z) + H 2 (z))/0.5(H l (z)-H 2 (z)) as shown in Figure 10, 

where an ordinary adaptive FIR filter may also be used to estimate the ratio of combined overall 
transfer functions. Its estimate may then be applied to the signal and noise period to cancel the noise 
by using estimated noise statistics. 



Figure 10. Analysis of identification of acoustic path transfer functions by adaptive FIR filter 
in beamforming structure during the periods of (a) noise alone and (b) signal plus noise. 




Based on this, an inverse kepstrum filter is applied in front of the adaptive filter as a rear-end 
application from the sum-and-subtract function from the beamforming structure as shown in 
Figure 11(a), where the inverse kepstrum filter estimates the polynomial denominator part, 
l/0.5(#i(z) -H 2 {zj) and the cascaded adaptive FIR filter estimates the polynomial numerator part, 
0.5(H l (z) + H 2 (z)) from the ratio of combined overall transfer functions. It is compared with the 
kepstrum method as shown in Figure 11(b), where the kepstrum filter estimates the minimum phase 
term only from the numerator polynomial part, 0.5(H ^z) + H 2 (z)) and the cascaded adaptive FIR filter 
estimates the remaining part from the numerator polynomial part, 0.5{H l {z) + H 2 {z)) , where this 
numerator polynomial part is to be an overall transfer function with a delay filter in the rear-end 
primary input, d n . 

Figure 11. Rear-end application of (a) inverse kepstrum filter and (b) kepstrum filter to 
beamforming structure. 
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To compare with the operation in an ANC structure, the same components of the acoustic transfer 
functions are used as Equation (29) and this is represented as an unknown system as in Equation (37). 
It is to be analyzed as shown in Figure 12, where the operation of the inverse kepstrum is to be 
compared with the kepstrum filter as a rear-end application to the sum-and-subtract function of the 
beamforming structure. 

From the unknown system, we have the numerator polynomial part from the ratio of overall 
combined transfer functions and it is represented as an overall transfer function and the denominator 
polynomial part works as a delay filter, described in Equation (37): 

H{z) = 0.5[(l + 2z _1 ) + (l + 0.2z _1 )] =l + l.lz _1 with one sample delay, D 1 (37) 



Figure 12. Block diagram for the operation comparison of an inverse kepstrum filter and a 
kepstrum filter as rear-end applications to a beamforming structure. 
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For the operation of inverse kepstrum filter K (z) and adaptive FIR filter W x {z) , each one 
is estimated as: 

K\z) = 1 and W x {z) = 1 + l.lz 1 with one sample delay, D 1 (38) 

This indicates that the rear-end inverse kepstrum k ' (z) works as a whitening filter with one sample 
delay and a cascaded adaptive FIR filter W x (z) estimated a numerator part from the ratio of combined 
overall transfer functions. On the other hand, for the operation of the kepstrum filter K(z) and the 
adaptive FIR filter W 2 (z) , each one is estimated as: 

K(z) = 1 + 0.9z _1 and W 2 (z) = 1 + 0.2z~ l with one sample delay, D~ [ (39) 

This indicates that the rear-end kepstrum estimates the minimum phase term only from the 
polynomial numerator part from the ratio of combined overall transfer functions, where the 
nonminimum phase term is reflected to the minimum phase term by the reciprocal polynomial as 
z~"H(z~ l ) . The cascaded adaptive FIR filter estimates the remaining term from the polynomial 

numerator of the ratio of combined overall transfer functions. Based on this, the unknown system 
Equation (37) may be estimated by the operations (38) and (39) as rear-end applications of the inverse 
kepstrum and kepstrum methods to the sum-and-subtract function of the beamforming structure, 
respectively. 
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Let us now check with the inverse of the overall transfer function, (33), where nonminimum phase 
term in denominator polynomial is no longer exist and overall transfer function is now obtained from 
numerator polynomial part of the ratio of overall combined transfer function. 

Inverse kepstrum and cascaded adaptive FIR filter estimates as: 

K\z) = 1 and W x {z) = 1 + l.lz 1 with one sample delay, D » 1 (40) 

where inverse kepstrum filter works as whitening filter and adaptive FIR filter estimates numerator 
polynomial part of the ratio of overall combined transfer function. 

On the other hand, kepstrum and cascaded adaptive FIR filter estimates as: 

K(z) = \ + 0.9z~ l and W 2 (z) = 1 + 0.2z~ l with one sample delay, D 1 (41) 

where kepstrum filter works as minimum phase filter and adaptive FIR filter estimate remaining term 
from the numerator polynomial part of the ratio of overall combined transfer function. 

It shows that the estimates, Equations (40) and (41) by both methods on the inverted transfer 
function are same as Equations (38) and (39), which are the estimates by the both methods on the 
direct transfer function. It indicates that both kepstrum and inverse kepstrum methods do provide a 
solution in beamforming structure, which may give a practical solution in a reverberant noise with 
nonminimum phase component from overall transfer function because it does not need to swap the two 
microphones position. The detailed analysis on nonminimum phase transfer function has been 
investigated between ANC and beamforming structures [27]. 

5. Experiments 

Experiments were implemented in both simulation tests on pc software and real tests using real 
nonstationary noise in a room environment. According to the main three considerations (signal 
distortion in the desired signal, noise reduction performance in noise with nonminimum phase 
components and convergence level in estimates of noise statistics with the use of a small amount of 
adaptive FIR filter weights), the performances achieved when using an inverse kepstrum filter were 
verified in both the ANC and beamforming structures. Furthermore, the rear-end application of the 
inverse kepstrum method in beamforming structure was also compared with two front-end applications 
of the kepstrum method [14] and the inverse kepstrum method [22] in beamforming structure. 

The methodology is based on the fact that the coefficients (kepstrum and inverse kepstrum) and 
weights (adaptive FIR filter) are continuously updated during noise alone period to estimate noise 
statistics. When the desired signal is applied to noise, the last updated coefficients and weights are 
frozen and applied to the desired signal and noise. For a precise test, it is programmed to stop updating 
coefficients and weights, and then these are applied to the desired signal and noise period. To check 
the strength in amplitude and distortion status in desired signal, simple three sine waveforms are added 
and used as the desired signal for both simulation and real tests. For the test using real noise, we use a 
nonstationary music sound, tuned to a certain radio station. For the use of an adaptive FIR filter, an 
NLMS algorithm has been used with the use of step size /u = 0.001 for the simulation test and ju = 0.5 
for the real test. For the processing, 2,048 frame size, sampling frequency of 22,050 Hz and Nyquist 
frequency of around 11,000 Hz have been chosen. Two preamplifiers and two microphones of 
unidirectional electret condenser type are used, placed 7 cm distance apart in broadside configuration 
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for the real test in a room [3.8 m(d) x 3 m(w) x 2.8 m(h)]. The performance is to be verified by 
comparing signal waveforms in time domain, spectra in frequency domain, estimated coefficients and 
weights arrays and its pole-zero placements. 

5.1. Simulation Test 

For the simulation test, the acoustic transfer functions of (29) are used as the unknown system, 
which has nonminimum phase component in noise. The desired signal, consisting of three frequencies 
(500 Hz, 550 Hz and 700 Hz), has arbitrarily been used as shown in Figure 14(a). 

5.1.1. Adaptive FIR filter in ANC and beamforming structures 

The first test is to verify the noise cancelling performance in the ANC structure by applying three 
adaptive FIR filter weights for the noise characteristic with nonminimum phase component in the 
polynomial numerator (29) and nonminimum phase component in the polynomial denominator (33) in 
the acoustic transfer function. It is also verified in the beamforming structure. 

From the simulation test based on the block diagram (Figure 9) of the ANC structure [Figure 7(a)], 
it is found that the noise spectrum with nonminimum phase term in the polynomial denominator (33) 
shows worse and much different performance than one of nonminimum phase term in the numerator 
(29), as shown in Figure 13(a). On the other hand, from the simulation test based on the block diagram 
(Figure 12) of the beamforming structure [Figure 10(a)], spectra of both Equations (29) and (33) show 
good and almost same spectral performance, as shown in Figure 13(b). From the test result, it is shown 
that the adaptive FIR filter works well with the use of a small amount of weights and shows 
consistency in spectra in the beamforming structure for the noise with both nonminimum phase cases, 
Equations (29) and (33). On the other hand, based on the test result in Figure 13(a), a fluctuation in 
spectrum is expected in ANC structure in the case that noise statistics is frequently changed in a 
reverberant environment. 



Figure 13. Comparison in spectra between direct transfer function (29) and its inverse 
transfer function (33) on application of adaptive FIR filter in (a) ANC structure and 
(b) beamforming structure. 




For the second test, the acoustic transfer function Equation (29) is estimated by using three adaptive 
filter weights and then it is applied to signal and noise. To get the best result, it is found that it is 
needed to set the three sample delay as D~ 3 to reduce signal distortion in the ANC structure. On the 
other hand, it is found that there is no need to set the delay in the beamforming structure. From the test 
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result, by using a small amount of three adaptive FIR filter weights, it is shown that the performance in 
the beamforming structure provides no signal distortion in the desired signal without any delay 
adjustment and also better noise reduction in noise with nonminimum phase term than the performance 
in the ANC structure, as shown in Figure 14. 

Figure 14. Simulation test based on 3 adaptive FIR filter weights: (a) waveform of original 
desired signal; (b) waveform of desired signal with three samples delayed in the ANC 
structure; (c) waveform of desired signal with no sample delay in the beamforming 
structure; (d) corresponding spectra of (b) and (c). 




From the above two test results in terms of noise reduction performance in reverberant noise with 
nonminimum phase (Figure 13) and signal distortion in the desired signal (Figure 14), it is found that 
an adaptive FIR filter works better in a beamforming structure than in an ANC structure. 

5.1.2. Application of inverse kepstrum method to two structures and comparison with kepstrum 

The objective is to verify the performance of the inverse kepstrum method in ANC and beamforming 
structures. For the test in the ANC structure, the acoustic transfer functions of (29) have been used, hence 
the unknown system is expected to be estimated as (30). As verified in Table 2 (a)-(iii) and Figure 15 (a), 
an adaptive FIR filter using one zero does not approximate to (30). It can be approximated by 
increasing the adaptive FIR filter weights to four from (a) so that it gives almost the same performance 
as Equation (30), as shown in Table 2 (b)-(iii) and Figure 15(b). On the other hand, by applying two 
inverse kepstrum coefficients to a cascaded adaptive FIR filter using one zero, it also gives almost the 
same performance as Equation (30) as shown in Table 2 (c)-(iii) and Figure 15(c). For the comparison, 
with the use of same size of two kepstrum coefficients, it indicates that four adaptive filter weights are 
needed to approximate (30) as shown in Table 2 (d)-(iii) and Figure 15(d). 

From the test result, it is analyzed that in the ANC structure, the application of an inverse kepstrum 
filter works well with an adaptive FIR filter in terms of convergence with a smaller amount in adaptive 
filter weights, rather than when the adaptive FIR filter is used with application of a kepstrum filter. 
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From the test results between (c)-(ii) and (d)-(ii) in Table 2, it indicates that inverse kepstrum method 
needs only 50% of the adaptive FIR filter weights needed by the kepstrum method. 

Table 2. Weights and coefficients arrays showing each estimate output from the simulation 
test based on the block diagram in Figure 9. 



(a) Adaptive filter weights (one zero) 
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1.712 












(b) Adaptive filter weights (three zeroes) 
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1.800 


-0.359 


0.069 
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(d) (i) 
(ii) 


Kepstrum coefficients (two poles) and 
Adaptive filter weights (three zeroes) 


(0 
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0.299 


0.044 










(ii) 
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1.501 


-0.846 


0.234 








(iii) 
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1.800 


-0.352 


0.048 


0.032 


0.010 





Note: iii=i*ii, where * indicates convolution. 



Figure 15. Snapshot of pole-zero placement: (a) one zero from adaptive FIR filter 
(b) increased to three zeros from (a) (c) two poles from inverse kepstrum coefficients and 
one zero from the adaptive FIR filter (d) two poles from kepstrum coefficients and increase 
to three zeros from the adaptive FIR filter. 
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For the test in beamforming structure, acoustic transfer functions of Equation (29) has been used, 
hence unknown system is expected to estimate as Equation (37). As verified in Table 3 (a)-(iii) and 
Figure 16(a), adaptive filter using one zero is well approximated to Equation (37). 



Table 3. Weights and coefficients arrays showing each estimate output from the simulation 
test based on the block diagram in Figure 12. 



(a) Adaptive filter weights (one zero) 




1 


1.101 
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(b) (i) Inverse kepstrum coefficients (two poles) and 






(ii) Adaptive filter weights (one zero) 








(i) 
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(d) (i) Kepstrum coefficients (two poles) and 








(ii) Adaptive filter weights (eight zeroes) 
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0.898 


0.403 










(ii) 


1 


0.200 
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-0.000 
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Note: (A) iii=i*ii, where * indicates convolution. 

(B) (d) (ii) and (iii) shows only seven weights estimates in arrays. 



Figure 16. Snap shot of pole-zero placement: (a) one zero from adaptive FIR filter (b) two 
poles from inverse kepstrum and one zero from adaptive FIR filter (c) two poles from 
kepstrum filter and one zero from adaptive FIR filter (d) two poles from kepstrum filter 
and increased to eight zeros from adaptive FIR filter. 
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Figure 16. Cont. 




(c) (d) 



Based on this, by using two inverse kepstrum coefficients, it can be verified that it can also be 
approximated to Equation (37) as shown in Table 3 (b)-(iii) and Figure 16(b). On the other hand, by 
applying two kepstrum coefficients to the cascaded adaptive FIR filter using one zero, it cannot be 
approximated to (37), as shown in Table 3 (c)-(iii) and Figure 16(c). For the comparison, with the use 
of same size of two kepstrum coefficients, it indicates that eight adaptive filter weights are needed to 
approximate (37), as shown in Table 3 (d)-(iii) and Figure 16(d). 

From the test result, it is also analyzed in a beamforming structure and shows that application of 
inverse kepstrum filter works well with the adaptive FIR filter in terms of convergence with much 
smaller amount in adaptive filter weights, rather than when the adaptive FIR filter is used with 
application of a kepstrum filter. From the test results between (b)-(ii) and (d)-(ii) in Table 3, it 
indicates that inverse kepstrum method needs only 25% (75% less) of the adaptive FIR filters than the 
kepstrum method. 

From the comparison results between inverse kepstrum filter and kepstrum filter to adaptive FIR 
filter in terms of convergence in noise statistics and its pole-zero placements, it is found that the 
application of the inverse kepstrum filter could give a convergence benefit with the use of a small 
amount in adaptive FIR filter weights in the ANC structure [Table 2 (c)-(ii)], and is much more 
effective in a beamforming structure [Table 3 (b)-(ii)]. 

5.2. Real Tests 

The simulation test results suggest that the inverse kepstrum should achieve more noise reduction 
without signal distortion in the desired signal by using small amount of adaptive FIR filter for the real 
tests in a realistic reverberant environment. Therefore, the inverse kepstrum has been tested in a 
beamforming structure [Figure 11 I-(a)] in a room. To verify the performance in signal, the desired 
signal [Figure 17 I-(a)] consisting of three frequencies (500 Hz, 550 Hz and 700 Hz) has been used, 
and music sounds tuned to a radio station [Figure 17 I-(b)] has been used as real nonstationary noise in 
a room. From the test on the ANC structure, it is shown that the shape of the signal waveform has been 
distorted in time domain, as shown in Figure 17 I-(c). On the other hand, it is shown that almost the 
same shape of the signal waveform has been maintained in the beamforming structure as shown in 
Figures 17 I-(d). The performance has also been compared in spectra in frequency domain as shown in 
Figure 17 II-(e). In the beamforming structure, the performance has been compared between 200 
adaptive FIR filter weights alone [Figures 17 II-(f)-(i)] and 32 inverse kepstrum coefficients with 50 
cascaded adaptive FIR filter weights [Figure 17 II-(f)-(ii)], where the test results show that inverse the 
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application of the kepstrum with reduced adaptive FIR filter weights gives even better noise reduction 
performance while maintaining the original desired signal. It indicates that the application of the 
inverse kepstrum could provide a benefit of more computational simplicity than a kepstrum in 
processing as well as the ordinary adaptive FIR filter (Table 1). To reduce the computational 
complexity, the alternative method could be considered by decreasing the window size but the problem 
has been found with the widening sidelobe on the desired signal. By using the inverse kepstrum, it has 
been found that increasing to 4096 the window size is acceptable for the real-time processing with the 
narrowing sidelobe on the desired signal [21]. 

Figure 17. Real test in a room environment: (I) — waveforms of (a) desired signal (b) real 
nonstationary noise (c) processing with three samples delayed in ANC structure 

(d) processing with no delay in beamforming structure, (II) — corresponding spectra from 

(e) I-(c) and I-(d) (f) comparison between (i) processing with 200 adaptive FIR filter 
weights alone and (ii) processing with 32 inverse kepstrum coefficients and 50 adaptive 
FIR filter weights in beamforming structure. 




H-(e) II-(f) 



Note: O in (c) and (d) indicates a distorted part of desired signal (a). 

Furthermore, the performance [Figure 17 II-(f)-(ii)] of the inverse kepstrum in rear-end 
beamforming [Figure 1 1(a)] has been compared with the front-end application of the inverse kepstrum 
[Figure 18(a)] and kepstrum [Figure 18(b) in a beamforming structure. With the use of a kepstrum (or 
inverse kepstrum) coefficients to 32 and a reduced size of adaptive filter weights to 50, it is now to 
verify the performance on signal distortion in a desired signal when the desired signal is applied to the 
noise estimates of coefficients and weights, which are abruptly frozen on the rapid change of noise 
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statistics. To verify the performance, the average discrepancy between strengths of three frequencies 
after processing has been measured since each input frequency is applied with equal strengths in dB. 
For a trial test, the same desired signal is applied at random intervals to the frequently changing real 
noise environment, where the estimate of noise transfer function is randomly frozen accordingly. 
Figure 18 shows snapshot of a performance among the application of front-end kepstrum, front-end 
inverse kepstrum with comparison of application of rear-end inverse kepstrum in time and frequency 
domains. 

Figure 18. Front-end application of: (a) inverse kepstrum and (b) kepstrum in 
beamforming structure; (c) waveform in the time domain by inverse kepstrum; 
(d) waveform in time domain by kepstrum; (e) spectra in frequency domain between 
Figure 17 II-(f)-(ii) and (c); (f) its spectra in frequency domain between Figure 17 II-(f)-(ii) 
and (d); [commonly based on 32 inverse kepstrum coefficients (or 32 kepstrum 
coefficients) and 50 adaptive FIR filter weights]. 
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O in (c) and (d) indicates a distorted part from desired signal (Figure 17 I-(a)) in time domain. 

i — i in (e) and (f) indicates that distorted part in desired signal from the application of rear-end inverse kepstrum 
(Figure 17 II-(f)-(ii)) is compared with one of application from (c) front-end inverse kepstrum and (d) front- 
end kepstrum in frequency domain. 



As shown in Figure 18(c,e), the test results show that front-end inverse kepstrum provides better 
noise reduction but with attenuated amplitude in a desired signal. On the other hand, the front-end 
kepstrum [14] shows more strength in signal amplitude with almost same noise reduction as shown in 
Figure 18(d,f). However, it has been found from the test that both front-end methods [Figure 18(a,b)] 
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are more vulnerable to signal distortion than the rear-end inverse kepstrum method in beamforming 
structure. As shown in Table 4, the level of signal distortion has been compared by measuring an average 
discrepancy of signal strength in dB from front-end kepstrum [14], front-end inverse kepstrum [22] and 
rear-end inverse kepstrum, where it has been calculated as 0.25 dB, 0.45 dB and 0.15 dB, respectively. 



Table 4. Comparison of average discrepancy on application of front-end kepstrum [14], 
front-end inverse kepstrum [22] and rear-end inverse kepstrum. 



Application 


Front-end kepstrum 


Front-end inverse 
kepstrum 


Rear-end inverse 
kepstrum 


Frequency (Hz) 


500 


550 


700 


D 


500 


550 


700 


D 


500 


550 


700 


D 


Trial 1 


30 


29.5 


29.5 


0.25 


20 


20.5 


21 


0.5 


25 


25 


24.5 


0.25 


Trial 2 


29.5 


30 


30 


0.25 


21 


20 


20 


0.5 


25 


25 


25 


0 


Trial 3 


29.5 


30 


29.5 


0.25 


20 


20 


20.5 


0.25 


25 


25 


25.5 


0.25 


Trial 4 


30 


29.5 


30 


0.25 


21 


20 


21 


0.5 


25 


25 


25 


0 


Trial 5 


29.5 


29.5 


30 


0.25 


21 


20.5 


20 


0.5 


25 


24.5 


25 


0.25 


Average D (dB) 


0.25 


0.45 


0.15 



Note: D indicates discrepancy in dB from signal strength of desired signal in frequency domain. 



5.3. Summary 

For real-time processing in a realistic reverberant environment, the kepstrum method and inverse 
kepstrum method have been applied to ANC and beamforming structures, hence its performance on 
modified application [4,6-8] to desired signal and adaptive filter has been investigated in detail in 
this paper. 

Firstly, the reverberant nature of most rooms gives rise to nonminimum phase components in the 
acoustic transfer function [28] and its inverse is often required in a realistic reverberant 
environment [26], hence it might produce instability in processing. Therefore, the acoustic transfer 
function and its inverted transfer function have been tested in ANC and beamforming structures. 
Secondly, the acoustic noise transfer function changes rapidly and frequently in a realistic reverberant 
environment and the estimated noise statistics for the acoustic transfer function are abruptly frozen. It 
is then instantly applied during the desired signal period, which might cause a signal distortion. To 
verify the performance, a discrepancy in dB has been compared among the application of front-end 
kepstrum, front-end inverse kepstrum and rear-end inverse kepstrum in a beamforming structure. 
Thirdly, a fast convergence in adaptive filter application is essential in real-time processing so that a 
fast convergence with small amount of adaptive filter weights has been tested in both kepstrum method 
and inverse kepstrum method. 

For an accurate discrepancy measurement in a desired signal, three added simple sine waveforms 
(disregarding whether it is narrowband or wideband in this test) have arbitrarily been used as a sample of 
a desired signal so that the distorted amount in a desired signal could be measured by calculating the 
average discrepancy in dB to check the consistency level of the signal strengths in dB from the original 
desired signal. For the precise testing in estimate of coefficients and weights, instead of using automatic 
VAD, it has been programmed to stop and make the last updated coefficients and weights to be frozen 
on demand of the application, and then it is applied to the desired signal and real noise period. 
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From the test results on main three considerations using the above mentioned methodology, it can 
be summarized that for the application of adaptive FIR filter to noise cancelling scheme, it is found 
that adaptive filter works better in a beamforming structure than the ANC structure in terms of signal 
distortion in desired signal, noise reduction in noise with nonminimum phase component (Figure 14) 
and consistency of estimate in noise statistics on its inverted transfer function (Figure 13). 

For the rear-end application of the inverse kepstrum to beamforming structure, the inverse kepstrum 
method gives better convergence with a much smaller amount of adaptive FIR filter weights than the 
kepstrum method (Table 3), which contributes to better real-time processing (Table 1). It has also been 
found that the combination of rear-end inverse kepstrum and cascaded adaptive FIR filter gives better 
noise reduction performance with highly reduced adaptive filter weights than the sole application of an 
adaptive FIR filter (Figure 17). In addition, it gives less signal distortion in the desired signal than 
the front-end applications of kepstrum method and inverse kepstrum method (Figure 18, Table 4). 

6. Conclusions 

An adaptive FIR filter has shown a performance distinction between ANC and beamforming 
structures in terms of signal distortion in the desired signal and nonminimum phase in noise, where the 
beamforming structure provides better performance than the ANC structure on application of an 
adaptive FIR filter. Based on this, the innovations-based inverse kepstrum method has been applied to 
the beamforming structure, and its performance has also been compared with the kepstrum method and 
inverse kepstrum method as front-end applications in the beamforming structure. The simulation and 
real tests show that the innovations-based whitening inverse kepstrum method has provided a 
promising result in rear-end beamforming structure in terms of signal distortion in the desired signal, 
noise reduction in noise with nonminimum phase and convergence level in estimation of noise 
statistics with the use of a small amount of adaptive FIR filter weights. Furthermore, the inverse 
kepstrum method provides the computational simplicity in processing so that it could be a benefit to a 
practical real-time adaptive noise cancelling scheme. Further analysis and investigation of the inverse 
kepstrum method will be performed using ECG (electrocardiogram) signals in biomedical signal 
processing for noise cancelling schemes. 
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